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Abstract 

O 
O 

o 

\ In this paper the Stewart-Lyth inverse problem is introduced. It consists 

Q ' of solving two nonlinear differential equations for the first slow-roll param- 

eter and finding the inflaton potential. The equations are derived from the 
Stewart-Lyth equations for the scalar and tensorial perturbations produced 
! during the inflationary period. The geometry of the phase planes transverse 

to the trajectories is analyzed, and conclusions about the possible behavior 
' for general solutions are drawn. 

l> 

. I. INTRODUCTION 

l> 

Q \ The analysis of observations of the cosmic microwave background radiation, recently 

O ■ reported by the experiments Boomerang and Maxima-1 confirmed that, at present- 

^ • day, inflation |^,^ still remains as the favourite model for the origin of structure in the 
Universe. The main reason is that the currently obtained observational data on structures 
2 ! in the Universe are interpreted in a natural way within the framework of the Gaussian 
^ \ adiabatic and nearly scale-invariant density perturbations that the usual models produce. 
^ ' On the other hand, it is one of the simplest paradigms within which rigorous theoretical 

predictions can be achieved. 
^ . However, the inflationary paradigm is quite a broad one, and there are several equally 

satisfactory implementations of the inflationary idea. The simplest scenario arises when 
the dynamics of inflation (both classical and quantum) are dominated by a single scalar 
field (inflaton) evolving in a nearly flat potential. Even in this case, there are a large 
number of candidates for the potential. For this scenario it is well established 0,^] that, to 
a good approximation, the scalar and tensor perturbations will take on a power-law form, 
with the tensor ones giving a subdominant (and almost negligible) contribution. A way 
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of discriminating between candidates for inflationary potential is the comparison between 
observed and predicted perturbations spectra. However, the present level of accuracy of the 
observations is well below the accuracy of the predictions for most of the models [0 , allowing 
a large number of them to remain as candidates. 

Another way of accomplishing the task of choosing a proper potential is to reconstruct 
it from the observed spectra , but again the low accuracy of measurements as well as the 
fact that the reconstruction is based on a Taylor expansion of the potential (with higher 
order derivatives given in dependence of the amount and quality of observational data), 
implies a large uncertainty in the outcomes of the standard reconstruction procedure, from 
the point of view of the potential uniqueness |^. Anticipating the near-future launch of 
satellites capable of measuring microwave background anisotropics to an accuracy of a few 
percent or better, across a wide range of angular scales ^ , an assessment of the accuracy of 
predictions of the anisotropics for given cosmological models has been stressed. Grivell and 
Liddle |^ have confirmed that for most models of inflation, the Stewart and Lyth analytic 
calculations |jlO[ should give extremely accurate predictions. 



Lyth and Stewart [|ri|] began with the precise calculation for power-law inflation |]T2[ and 
then they went on to use this exact result to analytically compute the next-order slow-roll 
correction to the standard formula |TD[ . At this level of approximation, the Stewart and Lyth 



equations for the spectral indices [10| can be rewritten as nonlinear differential equations in 



terms of the first slow-roll parameter e [13,14]. Hence, a third option for the determination 



of the inflationary potential, that we shall discuss in this paper, is to use observational 
information about density perturbations and gravitational waves spectra as input in these 
differential equations, to solve for e and to find the corresponding inflationary potential. We 
call this procedure the Stewart-Lyth inverse problem. 

In this work, after introducing the Stewart-Lyth inverse problem, we test the feasibility of 
the method by finding the proper potential for the well-known scenario of power-law inflation. 
Nevertheless, the involved system of differential equations composed by a second order and 
a first order equations, have proved to be very difficult to solve for general functional forms 
of the perturbations spectra. Taking this into account, in this paper we restrict our analysis 
to a qualitative study of the dynamics described by these equations. To draw conclusions 
about the behavior of the general solutions, the equations are studied using the spectral 
indices as parameters. For the corresponding reduced equations, we analyze the phase space 
for the second order equation which presents a singularity that strongly determines the flow 
geometry. Not cyclical orbits are found. The reduced first order equation is solved for any 
value of the tensorial spectral index and the possible behaviors are analyzed with all of the 
solutions being monotonic. Except in the case of null tensorial index, the singularity is also 
observed. In general, for solutions of the Stewart-Lyth inverse problem with smoothly and 
slowly changing spectral indices, periodic, quasiperiodic, or chaotic general solutions should 
not be expected. The theoretical analysis as well as observations suggest that power-law 
solutions are still valid. 

In the next section we briefly describe the theoretical frame for the Stewart-Lyth cal- 
culations and present their originally algebraic equations in terms of the spectral indices. 
In Sec. ^ we introduce the Stewart-Lyth inverse problem. Section |rV| is devoted to the 
qualitative analysis of the phase-spaces of the reduced second order equation for the first 
slow-roll parameter. The analysis of the solutions for the reduced first order equation is 
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presented in Sec. 0. We summarize the main results obtained in Sec. 



ivn 



II. THE STEWART-LYTH EQUATIONS 

A. The single scalar field scenario 

The theoretical frame for the Stewart-Lyth calculation is the flat Friedmann-Robertson- 
Walker universe containing a single scalar field equivalent to a perfect fluid with equations 
of motions given by 




+ 3H<j) = -y' (0) 



(1) 
(2) 



where (p is the inflaton, V{(f)) is the inflationary potential, H = a/ais the Hubble parameter, 
a is the scale factor, dot and prime stand for derivatives with respect to cosmic time and 
respectively, k = Svr/mpi is the Einstein constant and mpi is the Planck mass. 

In this framework, the first three slow-roll parameters were respectively defined in 
Ref. [1131, we shall write them here as in Ref. 
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Up to a constant, the first slow-roll parameter (|^) is a measure of the relative contribution 
of the kinetic energy to the total energy of the field. By definition e > 0, and because inflation 
can be defined as a > 0, for inflation to proceed: e < 1. 



B. Power-law inflation 

Few models of inflation allow to exactly calculate scalar and tensorial perturbations. One 
of them is the power-law model [T^, a particular scenario of inflation, where 



a{t) oc fP , (6) 

i/(0)cxexp(^-^0^, (7) 

^(0)«exp(^-^0^, (8) 

with p being a positive constant. It follows from Eqs. (^), (^), and (^) that in this case the 
slow-roll parameters are constant and equal each other. 
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C. The indices equations 



Exact expressions for the asymptotic scalar (density perturbations) and tensorial (gravi- 
tational waves) power spectra for the case of power-law inflation was correspondingly derived 



by Lyth and Stewart [|TT| and Abbot and Wise |I6 . 

Assuming that the deviation of the higher slow-roll parameters from e is small (power-law 
approximation) and that e is small with respect to unity (slow-roll approximation), Stewart 



and Lyth |TD[ derived next-to-leading order expressions for both spectra. In terms of the 



spectral indices, these expressions are 

1 - Usik) c^Ae-2ri + 8{C + l)e^ - (IOC + 6)eri + 2C^^ (9) 
ririk) ^ -2e [1 + (2C + 3)e - 2(C + l)r]] , (10) 

where the notation is that of [0]; ns{k) and rixik) are the scalar and tensorial spectral indices 
respectively, k is the wave number of the comoving scale, and C ~ —0.73 is a constant 
related to the Euler constant originated in the expansion of Gamma function. The symbol 
~ is used to indicate that these equations were obtained using the power-law and slow-roll 
approximations. Hereafter we shall use the = sign in our calculations, but the meaning of 
approximation should be added whenever it applies. 

For a giving expression of the scale factor, the Hubble parameter and the potential 
are determined, and then by substituting definitions (^, (H), and (||) in the Stewart-Lyth 
equations (|^) and (|10|), the scale-dependent spectral indices are obtained. For instance, 
giving Eq. (^ and substituting in Eqs. (^, (|), and (|]), one obtains € = = = 1/p, which 
in turn substituted in Eqs. (^ and (|T^) yield — 1 = rij- = const. This is an alternative 
definition for power-law inflation. We should note that the spectral indices are directly 
related with the power spectra (for details see Ref. M). 



III. FORMULATION OF THE STEWART-LYTH INVERSE PROBLEM 

In this section we shall formulate a method for finding the infiaton potential using ob- 
servational information on the spectral indices. Denoting T = 0^/2 and using definitions 
(0), (H), (^, together with Eqs. (|I]) and (||) then, in a straightforward manner, we obtain 



T kT 



T + v m 

dT 
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(11) 

(12) 



Defining r = In H^, 5{k) = nT{k)/2 and A(fc) = [n,(fc) - l]/2, and substituting Eqs. (0), 
(|l^), and (|T^) in Eqs. (H) and ([T0|) , the indices equations in terms of the first slow- roll 



parameter e and its derivatives with respect to r (e = de/dr and e = d'^e/dr^) become 

2Cee" - (2C + 3)ee - e + + e + A = , (14) 
2(C + l)ee- 6^-6-5 = 0. (15) 
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The second order nonhnear differential equation (|T4D was first introduced in Ref . [0 , while 



the first order equation (ITq) was derived first in Ref. [14|. Notice that Eqs. (14) and (15) 



are just alternative representations of the Stewart-Lyth Eqs. {^) and ([T0|). The approach 



of Ref. for the potential reconstruction is incomplete because is restricted to the case 
A{k) = const, and, on the other hand, the equation for the tensorial spectral index (|I5|) 
is considered only in the trivial case of power-law inflation [|14|. As we shall see, Eq. (plSf ) 
imposes rigorous constrains upon the set of solutions of Eq. (|Tip. 

By the Stewart-Lyth inverse problem we mean the problem consisting in solving Eqs. (p!4D 
and (|T3p for e, given expressions for the spectral indices, and finding the corresponding 
inflaton potential using the definitions of the first slow- roll parameter (|^) and ([Tl|). We 
prefer this denomination instead of "reconstruction" because we will use explicit functional 
forms of the spectra as input in the problem and correspondingly we will obtain functional 
forms for the potentials rather than pieces of information about them. 

As was already mentioned, having an expression for e(r) the corresponding potential as 
a function of r can be obtained from Eq. ([Tl|): 

r(r) = i[3-e(r)]exp(r). (16) 

On the other hand, taking into account that (p = f(f) and using Eq. (^, the scalar field 
as a function of r is given by 

<^M = -^ Z^ + '^o, (17) 



'2k 



e r 



where (pQ is an integration constant. 

Finally, the inflationary potential as parametric function of the inflaton can be given: 



rx- (18) 




The above expressions are similar to those used in Ref. ||13l but we would like to stress that 
the functional form for e in expressions ([T6|) and ([T7|) must be solution of both Eqs. ([1^ 
and dH). 

As a simple example of a resolution of the Stewart-Lyth inverse problem, let us ana- 
lyze the case of power-law inflation, where the spectral indices are constants. Substituting 
Eq. (p!5| ) in Eq. ( p^ with constant A and 6, an algebraic equation for e is obtained. Since 
the coefficients of this equation are scale-independent, all the solutions for this algebraic 
equation are just of the form: e = 1/p with p = const. Substituting this expression for e in 



Eq. (|T^ and after some algebra, we obtain for the Hubble parameter 

/7(0) = expf-./^(0-0o)V (19) 



We can see that, up to a constant, Eq. (|19]) is equivalent to (|^). Now, substituting 
e = 1/p and in Eq. (|T^), we obtain. 



y(0) = ^ exp -J- (0 - 0o) , (20) 
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in complete correspondence with Eq. (P). 

On the relevance of taking into account the first-order equation, we want to remark that, 
although Eq. ( [I^ ) has a large number of solutions for a constant scalar index, once Eq. ( p!5[ ) 
with S = const was used as a first integral of Eq. (p!4D, the unique remaining solution is just 
the potential given by Eq. (|20|). 

Even for a scale-independent scalar spectrum, Eq. (0), when it is taken with no regards 
to Eq. (|T3p, has proved to be very difficult to be analytically solved; solutions have been 
found only for a few fixed values of A (see Refs. ||T^JT7[| and references therein). This obstacle 
is even harder to overcome when an scale-dependent scalar spectrum is assumed; even for 
the first order Eq. (|1^) it is hard to find an explicit solution for e for any value of 5 = const, 
and even more difficult for 6 = 6[T{k)]. 

In the present work, we will focus in the dynamical aspect of the Stewart-Lyth inverse 
problem, i.e., in the dynamics of the first slow- roll parameter determined by Eqs. (|l^ and 
(|r^). The problem of integrating for the infiationary potential is left to be done in a near 
future. 



IV. PHASE SPACES ANALYSIS OF THE REDUCED SECOND ORDER 

EQUATION 

For an understanding of the dynamics behind Eqs. ([Mf ) and (p!5|) , the best approach seems 
to be a qualitative analysis of the corresponding phase-spaces. Not yet having an explicit 
expression for A(r) and S{t), we use the following approach in order to draw conclusions 
about the dynamics described by Eqs. (|T^ and (|I^): if we consider A(r) and 6{t) as the 
forcing element in these equations, then we can assume the dynamics to be characterized by 
one more dimension. The (e, e) planes [(r, e) for the first order equation] corresponding to 
the different values of this new coordinate are transverse to the trajectories given by these 
equations. Having the phase-portraits on the planes (e, e) for any value of A = const in 
Eq. ( p!^ and the solutions (r, e) for any 6 = const in Eq. (|15|), the geometry of the surfaces 
along which the real trajectories spread out could be outlined, assuming slow variation for 
A(r) and 6{t). 

Hereafter, we will refer to Eqs. (0) and (|l^) with constant A and 6 as the reduced 



equations of the Stewart-Lyth inverse problem, and the space for solutions depending on r 
as the extended phase-space. 

Solutions for the reduced first order equation (^) will be studied in the next section. 
Let us now proceed with the analysis of the dynamics given by Eq. ([M])- This equation can 
be rewritten as 



With the change of variables 



Xi = e . 
X2 = e. 



(22) 



we obtain the system 
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Xi 
X2 



F{xi,X2) 
G {Xi,X2) 



X2 



X2 . 



(23) 



For A 

the vector field, holds at every point {xi,X2 



const, the condition for existence of solution []T^, i.e., the continuity condition for 

eV = 1R^\ {(0,X2) , Vxa}. On the other hand. 



G{xi,X2^) - G{xi,X2 



(2). 



2\C\ 



{2G + 3) + — 

Xi 



^(1) 

X2 



^(2) 
X2 



(24) 



where the upper indices denote any two different values of X2- Then, the uniqueness condi- 
tion, i.e., Lipschitz condition, holds for the same set V. Therefore, unique solution for the 
equations system (p3D certainly exists at any point in V. In a similar fashion, the differen- 
tiability of solution with respect to initial conditions and parameters of the system is also 
satisfied. Hence, we can use the results of qualitative theory of dynamical systems in the 
plane for the study of system (|23|) [[l^] . 

The set of singular points of system (^) is the union of the set of fixed points {F = G = 0) 
and the X2 axis. We should note that around point Ps = (0, A) the field is well defined along 
directions for which X2 ^ A faster than xi ^ 0. This case will be analyzed in detail in 

Sec. rvn. 



A. Fixed points 

The fixed points of this system are the set {{xi, X2)} such that 

F {xi,X2) = , 
G (xi, X2) = . 

Therefore X2 = 0. From G (xi, X2) = and for xi 7^ 0, we have 

xl + xi + A = 0. (25) 
The roots of this quadratic polynomial are 



1 y/T^lA 
a± = -2± 2 ' 



where the symbol a± stands for the two different roots of Eq. pq). These roots are placed 
symmetrically with respect to Xi = —0.5. From now on we shall use the symbol a without 
lower indices while referring to any of these roots. 

According to the Grobman-Hartman theorem |^0|, the local behavior near the equilib- 
rium position depends on the linear part (if exists) of the vector field, in those cases where 
the eigenvalues of the Jacobian matrix of the vector field associated to the system have 
nonzero real part. Hence, in order to develop local analysis of the equilibrium positions, we 
calculate 
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Because the equilibrium positions are of type (a, 0) the Jacobian matrix in these points 
is written as 



J(a,0) 



1 

9i 92 



where 



91 

92 



A 

2C 
1 



(28) 

(29) 
(30) 



The characteristic polynomial, corresponding to the eingenvalue problem of Eq. 

p(X) = - g2X - 91 , 

with roots 

1 



IS 



A. 



2^ 2 



9l 



4^1 



(31) 



(32) 



We study the possible behavior near the equihbrium positions depending on the values of 
the polynomial roots (^). For a complete discussion, definitions and main results, see Refs. 
||18|^0| . In general, one is concerned about whether the eigenvalues are real or complex, and 
if their real parts are greater, equal or less than zero. 

The roots (|32|) depend on the value of discriminant D = g2 + 4(7i; they will be real if 
D > 0. Hence, the following inequality holds: 



D 



c 



2CJ 



C 2C2 



1 

- + 

a 



This inequality can be written as 

"(2C + 3)2 - 8C 
The polynomial in Eq. (|34D has roots 



a 



2(2C + 3)a + 8CA + l > 0. 



(33) 



(34) 



(2C + 3) ± 2^2C (1 - [{2C + 3)2 - 8C] A) 
(2C + 3)2-8C ■ 



(35) 



We shall see soon that the fixed point (q;_, 0) is always a saddle point. Hence, the current 
analysis has sense only for a+. Comparing the expressions for this root given by Eqs. (^Gf) 
and (p5D, solving for A and substituting the value of C ~ —0.73, we find out that roots ( |32D 
are real if and only if A < 0.124273. 

Now we are able to classify the fixed points. For A > 0.25, case (a), there are not fixed 
points and the flow is deformed from a laminar one only near the point ps (see Fig. |l|, case 
(a) and Sec. p.V C| ). Hereafter, in the figures we represent the flow direction by means of 
arrows. 
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FIG. 1. Reduced phase portraits for A = 0.3, case (a), and for A = 0.25, case (b). 



1.5 



For A = 0.25, case (b), one of the eigenvalues has zero real part and hence the stabihty of 
the unique and degenerated fixed point (—0.5, 0) cannot be determined by hnearization. The 
remaining eigenvalue has positive real part, then the local behavior transverse to the center 
manifold is controlled by the exponentially expanding flows in the local unstable manifold. 
According to Center Manifold theorem for flows the center manifold (it could be not 
unique) is tangent at (—0.5, 0) to the center subspace that, in this case, coincides with the 
Xi axis. On the other hand, the unstable manifold is tangent to the unstable subspace, i.e., 
the line X2 = 92^1 with g2, for this case, being (2C + 1)/2C ~ 0.315068 [Fig. ||, case (b)]. 
We recall that each subspace is spanned by the eigenvectors corresponding to eigenvalues 
(p^). Further, for A < 0.25 and A 7^ we have two nondegenerate fixed points over the Xi 
axis symmetrically distributed with respect to Xi = —0.5. 




FIG. 2. Reduced phase portrait for A = 0.227694, case (c), and for A = 0.124273, case (d). 
Let us proceed with the fixed points classification paying attention to the leftmost one. 
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FIG. 3. Reduced phase portrait for A = 0, case (e), and for A = —0.2, case (f). 

Obviously, the value of «_ for this point is always negative for any A and hence gi > 0. 
This way, the leftmost fixed point is always a saddle. The rate |A+| / |A_| is equal to 1 for 
A = -0.25 + 2(C + l)/(2C + 3) ~ 0.227694, value of A that corresponds to a+ = -0.64935. 
For A greater (less) than 0.227694 the trajectories converge to the saddle point faster (slower) 
than they diverge from it. 

For the rightmost fixed point the situation is richer. First, from Eq. ( pO]) one can see that 
(72 > for any A < 0.25 and A 7^ 0. For < A < 0.25 we have gi < 0, then in the interval 
0.124273 < A < 0.25, case (c), the rightmost fixed point is an unstable focus (Fig. As 
can be observed in the figure, one of the stable separatrices of the saddle point seems to be 
directly spiraling out from the focus. For < A < 0.124273, case (d), the rightmost fixed 
point is an unstable node (Fig. ||). In the case A = 0, case (e), (72 < and we have a single 
saddle located at (—1, 0) (see Fig. |^). For A < 0, case (f), the rightmost fixed point becomes 
a saddle (Fig. ^) located on the right half-plane. 



B. Closed orbits 

Note that the sign of the vector field divergence of system (p3D does not depend on A, 
thus the value of the parameter is of no importance while applying the Bendixson criterion 
TP| , which states that no closed path can be found in plane regions where the vector field 



divergence is null or not undergoes sign changes. This criterion, applied to system (p3D, 
divides the plane in four regions without divergence sign changes: 

Vi = {{xi,X2);xi > 0,Vx2}, with diverging vector field, 

V2 = {ixi,X2); —0.64935 < Xi < 0,Vx2}, with a converging vector field, 

V3 = {(—0.64935, X2); VX2}, with null divergence, and 

V4 = {{xi,X2);xi < —0.64935, VX2}, with positive divergence. 

It is worthy to note that, even if both half-planes (xi < 0, Xi > 0) are connected through 
point Ps, this connection should take place in one and only one direction. This way, no 
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closed orbit could be found in regions Vi and P3 or entirely lying in V2 or V4. Any closed 
trajectory should belong to U U V4. Then, three possible scenarios with closed path 
could take place. First, a homoclinic orbit (i.e., a saddle connection) enclosing the second 
fixed point (focus or node, we will refer to it as The second possibility is a limit cycle 
to which converges all of the trajectories starting from a+ and from which spiral out the 
orbits. Finally, it could be possible to find a homoclinic orbit with an embedded limit cycle 
around a+. Each of these scenarios has two possible realizations depending on whether both 
fixed points lie in V2 or the saddle lies in V4 and a+ lies in V2- The realizations for the 
more general scenario (a homoclinic orbit with an embedded limit cycle) are schematically 
represented in Fig. |[ 






X2 








(0,A) 


















Xi 



FIG. 4. Scheme of a homoclinic orbit with an embedded limit cycle. The dashed curve repre- 



sents the limit cycle and the dot-dashed line, the trajectory given by Eq. 36. The dotted line is 
= -0.64935. 



One conclude from the left scheme in Fig. ^that not closed trajectory will be found with 
both fixed points lying in V2 because sign changes in the derivative for the vector field along 
Xi, at each side of axis X2 = 0, are not allowed by the first equation of system (|23|). With 
regards to the realization plotted on the right part of Fig. ^, we shall see in Sec. [IV C| that 
for this range of Xi there is a trajectory (represented by the dot-dashed line on the figure) 
that can be fairly approximated by 



X2 



(36) 



then no closed orbit can neither exists here. 



C. Flow near point ps 

As we have already mentioned, a sharp feature of these phase spaces is the collisions of 
the trajectories with the border xi = 0. Our system possess here a vector field discontinuity 
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that strongly determines the behavior of the trajectories near this border. Such a piecewise 
continuous map is a more general system than those piecewise differentiable maps previously 
reported and analyzed in literature (see, for example, Refs. pl|j2^ ). 




A particularly interesting situation is exhibited by the trajectories near the point Ps 
(0,A). First of all, let us note that no statement can be done about the existence or 
uniqueness of solutions of system (^31) at this point. To analyze the behavior of the flow in 
the neighborhood of this point we rewrite the system as 

(37) 

2C) A - - 1 + [(3 + 2C)xi + 1] ^""^ ~ . (38) 

Xi J 

We shall look for a solution 0:2 = A + f{xi) for the phase curves in the neighborhood of 
the above mentioned point. Differentiating this expression with respect to Xi, taking into 
account that £2 = X\dx2ldx\ = X2df/dxi, substituting in Eq. (^) for the approximation 
Xi ^ 1 and f{xi) ^ 1 (x2 ~ A), we obtain the following nonhomogeneous linear equation 
for /(xi): 

^ _ J_L = (3 + 2C)A - 1 

dxi 2CAxi 2CA ' ^ ' 

For A 7^ 1/2C the solution of Eq. (|3|) is 



while for A = 1/2C, 



f{x,)={l-^^)x, + Kxr^\ (40) 



3 

f{xi) = —xi\\ixi + Kxi, (41) 

where K is an integration constant depending on the initial conditions. Hence, the behavior 
around the point Ps is described by the following family of curves 

^ ^ f A + (1 - ^) X, + KxT""^ , if A ^ 1/2C , .^2) 
' I 2^ + a^xilnxi + iTxi, ifA = l/2C. ^ ' 

Concerning the qualitative behavior of the trajectories near the point Ps there are four 
interesting intervals of A. First, for |A| ^ 1 the flow seems to ignore the existence of the 
special singular point (0, A). All the trajectories flows to (or from) the X2 axis along parallel 
lines with slope 1 + 3/2C —1.054794. From these lines only that intercepting the X2 
axis at ps could arrive to or depart from this point. For positive (but not too large) A the 
trajectories have the distribution observed in the left part of Fig. In this figure we present 
the case A = 1/(3 + 2C) ~ 0.64935 where the line with K = coincides with the stationary 
solution for Xi. Then, for A > one and only one trajectory could arrive to singular point 
Ps and one and only one trajectory could leave it. It is precisely this trajectory which can 
be represented by Eq. ( p6D for —0.64935 < a;+ < 0. For A = 0, in the left half-plane, all 
of the trajectories between the unstable separatrices of the saddle point and the ordinates 
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FIG. 5. Flow near the point ps for A > (A ~ 0.64935) and for A < (A ~ -0.2). 



0.002 „ 0.003 



axis move away from ps diverging from the line X2 = Xi. In the right half-plane, all of 
the trajectories converge asymptotically to the solution that behaves as X2 = Xi in the 
neighborhood of Ps [recall Fig. ^ case (e)]. All of the trajectories converge to Ps or diverge 
from it in the same direction. Finally, in the case of negative A (but again with not too 
large absolute value), the trajectories converge to ps or diverge from it tangent to the lines 
K = {) with slope given by Eq. p^) (Fig. ^, right part). 



V. THE REDUCED FIRST ORDER EQUATION 

We have already stressed the importance of Eq. ( [T3| ) for constraining the possible solu- 
tions of Eq. (p!4D. Let us consider now the equation 



2(C+l)ee-e2-e-5 = 



(43) 



with constant 5. 

The solutions of this equation also depend on the parameter value. For 5 > 0.25, we 
have 



In 



\e^ + e + 6\ 
B 



1 



: arctan 



e + 1/2 



^6 - 1/4 1/4/ C + 1 



where B is the integration constant. 

For 6 = 0.25, the solution of Eq. (|^) is that of the algebraic equation 



(44) 



[e + l/2)^-5exp 



C+1 e + 1/2^ 



0. 



(45) 



Now, for 6 < 0.25, and -0.5 - 0.5^1-45 < e(r) < -0.5 + 0.5^1-45, the solution is 
obtained from 
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„ /Vr^45-2e- 1\ ^™ 

— B , exp 

V VI -45 + 2e + 1 / 



C + 1 



(46) 



and, for 5 < 0.25 but e(r) < -0.5 - 0.5^1 - 45 and e(r) > -0.5 + 0.5^1 - 45, the solution 
results from 



+ e + 5 



^ /2e + 1 - Vr^45\ ^ 
— -D I , I exp 



r 



^2e + 1 + VI -45/ VC+1 
Note that the solution for 5 = is a special case of p7|), i.e., 



0. 



(47) 



e(r) = (eo + 1) exp 



r 



2{C+l) 



(48) 



where eg is the initial condition. 

Possible behaviors of solutions of Eq. ( ^Bj ) for different values of 5 are summarized in 
Fig. ^. The five interesting intervals are represented by vertical bands in the above plot. 




FIG. 6. Solutions of the reduced first order equation for different values of 5. In each of the five 
regions divided by vertical lines typical solutions for 5 in the corresponding interval are presented. 
In every region each curve branch that goes from left to right is a solution with different initial 
values. 



In each band the curves branch spreading from left to right is a solution starting from a 
different initial condition. We can see that all of the solutions are monotonic. In some cases 
they evolve unbounded, in some other cases bounded by the stationary solutions given by 
the roots of 
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VI. CONCLUSIONS 



We introduced the Stewart-Lyth inverse problem as the determination of the inflaton 
potential through expressions for the first slow-roll parameter obtained as a solution of 
differential equations. These equations were derived from the Stewart-Lyth equations for 
the spectral indices. We tested the feasibility of the method by solving the problem with 
constant spectral indices as input, corresponding with power-law inflation. 

To draw conclusions about the behavior of general r-depending solutions we analyze 
reduced equations using the spectral indices as parameters. The phase space for the reduced 
second order equation is richer than the generic ones, particularly due to the singularity 
e = and to the existence of the special singular point [0, 0.5{ns — 1)] near which the flow 
is deformed even when there are not fixed points. Do not exist cyclical orbits given by 
the reduced second order equation for any value of n^. The condition for the existence of 
stationary solutions with positive first slow-roll parameter is Ug < 1, and there is only one 
such a solution for every constant value of the scalar spectral index. 

The reduced first order equation was solved for any value of the tensorial spectral index. 
Five possible behaviors were found for the solutions: stationary, asymptotically increasing, 
asymptotically decreasing, boundlessly increasing, and boundlessly decreasing. The singu- 
larity e = was also observed, except in the case ut = 0. The condition for the existence of 
stationary solutions with positive first slow- roll parameter is < 0, and there is only one 
such a solution for every constant value of the tensorial spectral index. 

In general, for solutions of the Stewart-Lyth inverse problem with smoothly and slowly 
changing spectral indices in the expected range of values, we shall find that the trajectories 
would be confined in one of the sectors of the extended phase space divided by the axis 
e = 0. The exception could be those system with the tensorial index crossing through value 
ut = while the scalar index is greater than 1. Due to the lack of periodic solutions for 
the reduced equations, periodic, quasiperiodic, or chaotic extended solutions should not be 
expected. Suitable (and unique in each case) power-law solutions will exist if and only if 
Us < 1 and (consistently with the definition of power-law inflation) < 0. Taking this 
last result into account, estimations of the scalar index based in the recent observations 
by the collaborations Boomerang and Maxima-1 0,0] indicate that, despite its simplicity, 
power-law scenario is still a good candidate for the inflationary stage of the early universe. 

A great inside about general properties of scale-dependent solutions for smoothly and 
slowly changing spectral indices was gained in this work. Further efforts will be focus on 
obtaining solutions for the Stewart-Lyth inverse problem, i.e., explicit functional forms of 
the inflationary potential. 
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